Model atmospheres of irradiated exoplanets: The influence of stellar parameters, metallicity, and the C/O ratio

Many parameters constraining the spectral appearance of exoplanets are still poorly understood. We therefore study the properties of irradiated exoplanet atmospheres over a wide parameter range including metallicity, C/O ratio and host spectral type. We calculate a grid of 1-d radiative-convective atmospheres and emission spectra. We perform the calculations with our new Pressure-Temperature Iterator and Spectral Emission Calculator for Planetary Atmospheres (PETIT) code, assuming chemical equilibrium. The atmospheric structures and spectra are made available online. We find that atmospheres of planets with C/O ratios $\sim$ 1 and $T_{\rm eff}$ $\gtrsim$ 1500 K can exhibit inversions due to heating by the alkalis because the main coolants CH$_4$, H$_2$O and HCN are depleted. Therefore, temperature inversions possibly occur without the presence of additional absorbers like TiO and VO. At low temperatures we find that the pressure level of the photosphere strongly influences whether the atmospheric opacity is dominated by either water (for low C/O) or methane (for high C/O), or both (regardless of the C/O). For hot, carbon-rich objects this pressure level governs whether the atmosphere is dominated by methane or HCN. Further we find that host stars of late spectral type lead to planetary atmospheres which have shallower, more isothermal temperature profiles. In agreement with prior work we find that for planets with $T_{\rm eff}$ $<$ 1750 K the transition between water or methane dominated spectra occurs at C/O $\sim$ 0.7, instead of $\sim$ 1, because condensation preferentially removes oxygen.Comment: 30 pages, 20 figures. Accepted for publication in Ap

Planetary population synthesis coupled with atmospheric escape: a statistical view of evaporation

We apply hydrodynamic evaporation models to different synthetic planet populations that were obtained from a planet formation code based on a core-accretion paradigm. We investigated the evolution of the planet populations using several evaporation models, which are distinguished by the driving force of the escape flow (X-ray or EUV), the heating efficiency in energy-limited evaporation regimes, or both. Although the mass distribution of the planet populations is barely affected by evaporation, the radius distribution clearly shows a break at approximately 2 $R_{\oplus}$. We find that evaporation can lead to a bimodal distribution of planetary sizes (Owen & Wu 2013) and to an "evaporation valley" running diagonally downwards in the orbital distance - planetary radius plane, separating bare cores from low-mass planet that have kept some primordial H/He. Furthermore, this bimodal distribution is related to the initial characteristics of the planetary populations because low-mass planetary cores can only accrete small primordial H/He envelopes and their envelope masses are proportional to their core masses. We also find that the population-wide effect of evaporation is not sensitive to the heating efficiency of energy-limited description. However, in two extreme cases, namely without evaporation or with a 100\% heating efficiency in an evaporation model, the final size distributions show significant differences; these two scenarios can be ruled out from the size distribution of $Kepler$ candidates.Comment: Accepted for publication in ApJ; 24 pages, 16 figure

Observed spectral energy distribution of the thermal emission from the dayside of WASP-46b

We aim to construct a spectral energy distribution (SED) for the emission from the dayside atmosphere of the hot Jupiter WASP-46b and to investigate its energy budget. We observed a secondary eclipse of WASP-46b simultaneously in the g'r'i'z'JHK bands using the GROND instrument on the MPG/ESO 2.2m telescope. Eclipse depths of the acquired light curves were derived to infer the brightness temperatures at multibands that cover the SED peak. We report the first detection of the thermal emission from the dayside of WASP-46b in the K band at 4.2-sigma level and tentative detections in the H (2.5-sigma) and J (2.3-sigma) bands, with flux ratios of 0.253 +0.063/-0.060%, 0.194 +/- 0.078%, and 0.129 +/- 0.055%, respectively. The derived brightness temperatures (2306 +177/-187K, 2462 +245/-302K, and 2453 +198/-258K, respectively) are consistent with an isothermal temperature profile of 2386K, which is significantly higher than the dayside-averaged equilibrium temperature, indicative of very poor heat redistribution efficiency. We also investigate the tentative detections in the g'r'i' bands and the 3-sigma upper limit in the z' band, which might indicate the existence of reflective clouds if these tentative detections do not arise from systematics.Comment: 8 pages, 7 figures, accepted for publication in A&

On the massive young stellar object AFGL4176: High-spatial-resolution multi-wavelength observations and modeling

Deeply embedded and at distances of several kiloparsecs, massive young stellar objects (MYSOs) present numerous challenges for observation and study. In this work, we present spatially-resolved observations of one MYSO, AFGL 4176, together with survey and literature data, ranging from interferometric observations with VLTI/MIDI in the mid-infrared, to single-dish Herschel measurements in the far-infrared, and sub-millimeter data from APEX. We consider this spatially-resolved, multi-wavelength data set in terms of both radiative transfer and geometric models. We find that the observations are well described by one-dimensional models overall, but there are also substantial deviations from spherical symmetry at scales of tens to hundreds of astronomical units, which are revealed by the mid-infrared interferometric measurements. We use a multiple-component, geometric modeling approach to explain the mid-infrared emission on scales of tens to hundreds of astronomical units, and find the MIDI measurements are well described by a model consisting of a one-dimensional Gaussian halo and an inclined (\theta=60 deg) circumstellar disk extending out to several hundred astronomical units along a position angle of 160 deg. Finally, we compare our results both with previous models of this source, and with those of other MYSOs, and discuss the present situation with mid-infrared interferometric observations of massive stars.Comment: 15 pages, 14 figures. Accepted to Astronomy and Astrophysic

Observing transiting planets with JWST -- Prime targets and their synthetic spectral observations

The James Webb Space Telescope will enable astronomers to obtain exoplanet spectra of unprecedented precision. Especially the MIRI instrument may shed light on the nature of the cloud particles obscuring planetary transmission spectra in the optical and near-infrared. We provide self-consistent atmospheric models and synthetic JWST observations for prime exoplanet targets in order to identify spectral regions of interest and estimate the number of transits needed to distinguish between model setups. We select targets which span a wide range in planetary temperature and surface gravity, ranging from super-Earths to giant planets, and have a high expected SNR. For all targets we vary the enrichment, C/O ratio, presence of optical absorbers (TiO/VO) and cloud treatment. We calculate atmospheric structures and emission and transmission spectra for all targets and use a radiometric model to obtain simulated observations. We analyze JWST's ability to distinguish between various scenarios. We find that in very cloudy planets such as GJ 1214b less than 10 transits with NIRSpec may be enough to reveal molecular features. Further, the presence of small silicate grains in atmospheres of hot Jupiters may be detectable with a single JWST MIRI transit. For a more detailed characterization of such particles less than 10 transits are necessary. Finally, we find that some of the hottest hot Jupiters are well fitted by models which neglect the redistribution of the insolation and harbor inversions, and that 1-4 eclipse measurements with NIRSpec are needed to distinguish between the inversion models. Wet thus demonstrate the capabilities of JWST for solving some of the most intriguing puzzles in current exoplanet atmospheric research. Further, by publishing all models calculated for this study we enable the community to carry out similar or retrieval analyses for all planets included in our target list.Comment: 24 pages, 7 figures, accepted for publication in A&

Ground-based detection of the near-infrared emission from the dayside of WASP-5b

(Abridged) WASP-5b is a highly irradiated dense hot Jupiter orbiting a G4V star every 1.6 days. We observed two secondary eclipses of WASP-5b in the J, H and K bands simultaneously. Thermal emission of WASP-5b is detected in the J and K bands. The retrieved planet-to-star flux ratios in the J and K bands are 0.168 +0.050/-0.052% and 0.269+/-0.062%, corresponding to brightness temperatures of 2996 +212/-261K and 2890 +246/-269K, respectively. No thermal emission is detected in the H band, with a 3-sigma upper limit of 0.166%, corresponding to a maximum temperature of 2779K. On the whole, our J, H, K results can be explained by a roughly isothermal temperature profile of ~2700K in the deep layers of the planetary dayside atmosphere that are probed at these wavelengths. Together with Spitzer observations, which probe higher layers that are found to be at ~1900K, a temperature inversion is ruled out in the range of pressures probed by the combined data set. While an oxygen-rich model is unable to explain all the data, a carbon-rich model provides a reasonable fit but violates energy balance.Comment: 13 pages, 9 figures, accepted for publication in A&

Broad-band transmission spectrum and K-band thermal emission of WASP-43b as observed from the ground

(Abridged) We observed one transit and one occultation of the hot Jupiter WASP-43b simultaneously in the g'r'i'z'JHK bands using the GROND instrument on the MPG/ESO 2.2-meter telescope. From the transit event, we have independently derived WASP-43's system parameters with high precision, and improved the period to be 0.81347437(13) days. No significant variation in transit depths is detected, with the largest deviations coming from the i', H, and K bands. Given the observational uncertainties, the broad-band transmission spectrum can be explained by either a flat featureless straight line that indicates thick clouds, synthetic spectra with absorption signatures of atomic Na/K or molecular TiO/VO that indicate cloud-free atmosphere, or a Rayleigh scattering profile that indicates high-altitude hazes. From the occultation event, we have detected planetary dayside thermal emission in the K-band with a flux ratio of 0.197 +/- 0.042%, which confirms previous detections obtained in the 2.09 micron narrow band and Ks-band. The K-band brightness temperature 1878 +108/-116 K favors an atmosphere with poor day- to night-side heat redistribution. We also have a marginal detection in the i'-band (0.037 +0.023/-0.021%), which is either a false positive, a signature of non-blackbody radiation at this wavelength, or an indication of reflective hazes at high altitude.Comment: 14 pages, 9 figures, accepted for publication in A&

Spectral Evolution and Radial Dust Transport in the Prototype Young Eruptive System EX Lup

EX Lup is the prototype of a class of pre-main sequence eruptive stars defined by their repetitive outbursts lasting several months. In 2008 January-September EX Lup underwent its historically largest outburst, brightening by about 4 magnitudes in visual light. In previous studies we discovered on-going silicate crystal formation in the inner disk during the outburst, but also noticed that the measured crystallinity fraction started decreasing after the source returned to the quiescent phase. Here we present new observations of the 10 $\mu$m silicate feature, obtained with the MIDI and VISIR instruments at Paranal Observatory. The observations demonstrate that within five years practically all crystalline forsterite disappeared from the surface of the inner disk. We reconstruct this process by presenting a series of parametric axisymmetric radiative transfer models of an expanding dust cloud that transports the crystals from the terrestrial zone to outer disk regions where comets are supposed to form. Possibly the early Sun also experienced similar flare-ups, and the forming planetesimals might have incorporated crystalline silicate material produced by such outbursts. Finally, we discuss how far the location of the dust cloud could be constrained by future JWST observations.Comment: 12 pages, 4 figures, accepted for publication in Ap